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Abstract 

We have examined trapping of two-armed nearly vertical oscillations in 
polytropic disks. Two-armed nearly vertical oscillations are interesting in 
the sense that they are trapped in an inner region of disks with proper fre- 
quencies, if the inner edge of disks is a boundary that reflects oscillations. 
The frequencies of the trapped oscillations cover the frequency range of 
kHz QPOs to low frequency QPOs in LMXBs, depending on the modes 
of oscillations. Low frequency trapped oscillations are particularly inter- 
esting since their trapped region is wide. These low frequency oscillations 
are, however, present only when r(= 1 + 1/A'') is close to but smaller than 
4/3 (when spin parameter o* is zero), where is the polytropic index. 
The above critical value 4/3 slightly increases as a* increases. 

1 Introduction 

One of the possible candidates of high-frequency quasi-periodic oscillations (HFQ- 
POs) observed in low-mass X-ray binaries is resonantly excited p- and/or g-mode 
oscillations in deformed disks (Kato 2004, 2008a, b; Ferreira and Ogilvie 2008; 
Oktariani et al. 2010). The disk deformation required in this model is a warp 
or an eccentric deformation in equatorial plane. Recent numerical MHD simu- 
lations (Henisey et al. 2009) seem to have preliminary confirmed the presence 
of the excitation mechanism. 

The high frequency QPOs observed in black-hole candidates appear in pairs 
with frequency ratio of 3 : 2. This characteristic of high-frequency QPOs in 
black-hole candidates seems to be described by the above model, if we assume 
that geometrically thin disks are surrounded by hot tori and that the QPOs pho- 
tons are disk photons Comptonized in the tori (Kato and Fukue 2006) . The twin 
kHz QPOs observed in neutron stars, on the other hand, have always not the 3 
: 2 frequency ratio. They change their frequencies with time with correlation. 
If we want to describe such characteristics of kHz QPOs by the above model, a 
warp or eccentric disk deformation must have a time-dependent precession. In 
the case of neutron stars, distinct from the case of black holes, such precession 
of the deformation might be expected, since the central stars have surfaces and 
this might become causes of time-dependent precession through magnetic and 
radiative couplings between the central sources and the disks. It is not clear, 
however, whether possible time-dependent precession of deformation has a time 
scale consistent with the time variations of the kHz QPOs. 
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In this context, it is worthwhile examining whether there are other kinds of 
disk deformations that can become a possible source of resonant excitation of 
high-frequency oscillations in disks. This problem has been examined by Kato 
(2009), and it was suggested that two-armed nearly vertical oscillation^ can 
excite high-frequency p- and/or g-mode oscillations. Hence, it will be interesting 
to examine whether there are two-armed nearly vertical oscillation modes that 
can become a global deformation of disks. 

More interestingly, the two-armed vertical disk oscillations themselves may 
be origins of the varous types of QPOs observed in neutron-star LMXBs. This 
is because two-armed nearly vertical osillations occur in the inner region of disks 
and cover a wide range of frequency by difference of modes. 

Based on these considerations, we examine in this paper basic properties of 
the two-armed vertical disk oscillations. That is, we examine their characteris- 
tics of trapping, eigen-frequency and its dependence on polytropic index spec- 
ifying the vertical disk structure, using the mathematical formulations already 
prepared by Silbergleit et al. (2001) to study the trapping of the corrugation 
waves (c-mode oscillations). 

2 Overview of Trapping of Two- Armed Vertical 
Oscillations in Polytropic Disks 

Before having a rough image of nearly vertical oscillations in vertically poly- 
tropic disks, we briefly mention nearly vertical oscillations in isothermal disks, 
because it is instructive to have a rough image of wave trapping in disks. 

If a disk oscillation with azimuthal wavenumber m occur isothermally and 
locally in vertical direction in isothermal disks, the frequency, w, is given by 
(see, e.g., Okazaki ct al. 1987) 

(w - rnQf - 7in\ = 0, (1) 

where n denotes the node number (n = 1, 2, ...) of the oscillations in the vertical 
direction, and 57 and fl± are, respectively, the angular velocity of disk rotation 
at the radius in consideration and the vertical epicyclic frequency at the ra- 
dius. The oscillations are, however, cannot be purely vertical by inhomogeneity 
of the disk in the radial direction. They essentially have horizontal velocity 
components. In other words, the vertical oscillations described by equation ([1]) 
arc not rigorous. The vertical oscillations couple with purely horizontal inertial 
oscillations, (a; — mfl)^ — = 0, through pressure. If the coupling is considered 
with local approximations in the radial direction, we have a dispersion relation, 
which is written as (Okazaki et al. 1987) 

[{oj - mfl)^ - k'^][{uj - mn)^ - nfl^] = clk^iuj - m^^)^ (2) 

^ There will be no commonly accepted classification and terminology concerning disk os- 
cillations. Here, we classify oscillations into four types, i.e., p-mode, g-mode, c-mode and 
vertical p-mode oscillations (see, for example, Kato 2001; Kato et al. 2008). The oscillations 
considered here are the vertical p-mode oscillations. 
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where Cg is the acoustic speed, k, is the horizontal epicychc frequency, and k is 
the horizontal wavenumber of the oscillations. 

Since we are now interested in nearly vertical oscillations, i.e., (w — mfi)^ — 
nVl\ ~ (n = 1, 2, 3,...), equation ^ shows that their propagation region in 
the radial direction is described by 

(w - mVlf - nVL\ > 0, (3) 

since is always larger than k. 

We now proceed to polytropic disks with barotropic gas with polytropic 
index TV, i.e., p cx pt^+^Z^), where p and p are pressure and density, respectively. 
The frequency, w, of purely vertical oscillations is then given by (Silbergleit et 
al 2001, Kato 2005) 

(w - niQ)^ - «'„f7i = 0, (4) 

where is 

{1 for n = 1 

2 + 1/N for n = 2 (5) 
3 + 3/N for n^3. 

The case of isothermal disks corresponds to N = oo. 

In oscillations of n = 1, w — mQ is independent of N , but in oscillations of 
n ^ 1 it depends on A^. That is , for two-armed oscillations (m = 2), equation 
^ has solution given by 

w = 2f7 - (1 + r)i/2f^_L for n = 2 (6) 

and 

uj = 2n- (3r)i/2r2_L for n = 3, (7) 

where F = 1 + l/N, and F is usually in the range of F = 1 ^ 5/3. An 
interesting result shown in equation ([7]) is that oscillations with n — 3 can have 
low- frequencies for reasonable values of F. 

As in the case of isothermal disks, the nearly vertical oscillations are also not 
localized at a particular radius. They propagate in the radial direction. The 
propagation region is specified, as in the case of isothermal disks, by (see section 
4 for details) 

- mnf - «'„r2i > 0. (8) 

The propagation diagram showing the propagation region of the oscillations in 
the radial direction is given in figure 1 for oscillations with n — 1 and 2, and in 
figure 2 for oscillations with ti = 3. As the value of n increases, the propagation 
region of oscillations becomes narrower in the frequency - radius plane (i.e., 
propagation diagram), and finally for oscillations with ti = 3, the propagation 
region in the inner region of disks disappears for oscillations with F > 4/3. 

The next problem is to examine where the oscillations are really trapped 
and how much the frequency of the trapped oscillations is. To examine this 
problem, a careful treatment of the wave equation is necessary and this is done 
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in the subsequent sections. Here, however, in order to have a rough image of 
trapping, we discuss qualitatively what determines the frequency of the trapped 
oscillations. 

Let us assume tentatively that the frequency of the trapped oscillation is 
w and the outer boundary of the trapped region is r^. The inner boundary 
of the trapped region, ri, is taken at the inner edge of the disks, which is 3rg 
when the central source has no spin, being the Schwarzschild radius. Then, 
roughly speaking, Tc — ri must be on the order of the wavelength of the trapped 
oscillations. If some kinds of mean values of (a; — mfi)^ — and Cg in the 
trapped region are denoted, respectively, {{uj — m£l)'^ ~ k^) and (Cg), then the 
trapping condition is ((w — mfJ)^ — k^) ~ {Cs)/{rc — n)^ (for the fundamental 
mode in the radial direction), since oscillations propagate in the radial direction 
as inertial acoustic waves. This trapping condition is realized only when w has 
a particular value, say Wt- This is because for w > Wt, ((w — mfi)^ — K^)(rc — n)^ 
is too small compared with (Cg), since the width of the propagation region, i.e., 
Tc — n, is too narrow (see figures 1 and 2). On the other hand, for a; < Wt, the 
value of ((w — rnfi)^ — /t^)(rc — n)^ becomes larger than (Cg), since the width, 
Tc — Ti, becomes too large. To know the detailed value of Wt for the trapping, 
we must make numerical calculations, as is done in the subsequent sections. 

At this stage, without detailed numerical calculations, we can say that when 
r is large, the frequency of the trapped oscillations is low, if other parameters, 
including the radial distribution of Cso in disks, are fixed (see figure 4). Let us 
assume that the trapping condition has satisfied for a frequency, say Wt, when 
r has a value, say Ft. If F increases from Ft, rc and thus Vc — ri decreases 
from the value in the above case, if the frequency uj is fixed at ujt [see equations 
([6]) and ([7]) and curves in figures 1 and 2]. Thus, the trapping condition is not 
satisfied at w = Wt. For the trapping condition to be satisfied we must increase 
the width of the trapping region. This can be done by decreasing uj from LUt- 

3 Unperturbed Disks and Equations Describing 
Disk Oscillations 

3.1 Basic Assumptions and Unperturbed Disks 

We consider geometrically thin, relativistic disks. For mathematical simplicity, 
however, the effects of general relativity are taken into account only when we 
consider radial distributions of ri(r), ^(r), and ri±{r). They are, in turn, the 
angular velocity of disk rotation, the epicyclic frequencies in the radial and 
vertical directions. Except for them, the Newtonial formulations are adopted. 
Since geometrically thin disks are considered, fi, is approximated to be the 
relativistic Keplerian angular velocity, J7k(''), when its numerical values are 
necessary. Here, r is the radial coordinate of cylindrical ones [r^Lp,z), where 
the z-axis is perpendicular to the disk plane and its origin is the disk center. 
Functional forms of Q,]<i{r), K{r), and fl_L{r) are given in many literature sources 
(e.g., Kato et al. 2008) 
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The disks are assumed to consist of a barotropic gas with poly tropic index A'^, 
i.e., the pressure p and density p are related by p cx The hydrostatic 

balance in the vertical direction then gives (e.g., Kato et al. 2008) 



N 



po{r,z) = pQo{r)[l- , (9) 



,2 X 



Po(r,^)-Poo(r)(^l--^j , (10) 

where subscript represents the quantities in the equilibrium state and 00 are 
those on the equatorial plane. The acoustic speed defined by Cg = Tdpo/dpo is 
also given by 

cUr,z)^cUr)(l-j^y (11) 

where T is related to iV by T = (TV + 1)/7V or iV = l/(r - 1). The hydrostatic 
balance in the vertical direction also shows that the half-thickness H of disks is 
related to Cso and n± by 

n^H^ = 2Nc%. (12) 



3.2 Equations Describing Disk Oscillations 

We consider small-amplitude perturbations in the above disk. The perturba- 
tions are assumed to be proportional to exp[i{iot — mip)], where w and m are 
frequency and azimuthal wavenumber, respectively, of the perturbations. Then, 
hydrodynamical equations describing them are (e.g., Kato et al. 2008) 



iujpi + —^{rpQUr) - i—pou^ + -r-(poWz) = 0, (13) 
ror r oz 



dh\ 

iuJUr — 2Q.Uu, = 7^. (14) 

or 



dhi 
Oz 

where hi and uj are defined, respectively, by 



(16) 



/ii = — =Cs — , (17) 
Po P 

bj = UJ — mfl. (18) 

Here, {ur,u^,Uz), Pi, and pi are the Eulerian velocity, pressure, and density 
perturbations over the unperturbed ones, respectively. Equation (jl3p is the 
equation of continuity, equations (|14p - (jl6l) are in turn the r-, ip-, and z- 
components of equation of motion. 
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Elimination of Ur, u^p, u^, and pi from equations - PT]) leads to a partial 
differential equation of hi. In the processes, we neglect the radial variations of 
the unperturbed density and some other quantities, assuming that their charac- 
teristic radial scales are longer than the characteristic radial wavelength of the 
perturbations. Then, we have an equation describing the variation of hi in the 
form: 



Lj dr \ 



f- 



dhi 



dr 



dhi 



dz 



1 



hi = 0. 



(19) 



1 I d d 

ui^ Po dpa dz ' 

As is clear from the derivation processes, this equation is obtained under the 
framework of the Newtonian hydrodynamics. A fully general relativistic one 
corresponding to equation ([T9l) has been derived by Perez et al. (1997) and 
used by Silbcrgicit et al. (2001, 2008). In their equation the factor w inside 
and the factor l/w outside the first r-derivative d/dr in the first term on the 
right-hand side of equation p9|) are absent. However, relativistic factors related 
to metric appear in the first and third terms of equation (|19|) . These factors 
are on the order of unity, but not unity. For example, on the first term of 
equation ^9]). a metric factor g^^{r) is multiplied, which is 2/3 at 3rg (rg being 
the Schwarzschild radius) in the case of the Schwarschild metric. 

In this paper we do not take into account the relativistic effects except the 
functional forms of f2, k and u}± as mentioned before and start from equation 
(fT9|) . In the case of the polytropic gas mentioned in subsection 3.1, by changing 
the independent variables from (r, z) to (r, 77), where rj = z/H, we can reduce 
equation (fT9|) to 



or 



dhi 



dr 



1 



d 



(1 - r/2)JV 5^ 



1 



(20) 

The assumption of strong variation of perturbations in the radial direction 
ensures the WKB separability of variable in the above partial differential equa- 
tion in the form of hi{r,r]) = hi{r)g{r,ri). Here, the function g varies slowly 
with r. Extensive studies of disk oscillations by the WKB procedures have been 
made by Wagoner and his collaborators, e.g., Nowak and Wagoner (1992), Perez 
et al. (1997), Silbergleit et al. (2001). Here, we follow the procedures of Silber- 
gleit et al. (2001) used to study trapped c-mode oscihations, although notations 
adopted here somewhat different from theirs. 

By using the decomposition of hi (r, z) mentioned above, we can reduce equa- 
tion (|20|) into the following set of two ordinary differential equations (Silbergleit 
et al. 2001): 



1 



(1 - 772)^ dr] 



1 



2NKg = 0, 



and 



dr 



u> 



dhi 



dr 



2NKhi = 0, 



(21) 



(22) 



where 2NK is the separation constant. Since the term uj'^H'^ /c^q, which is equal 
to 2Noj'^ /il'^, in equation (|2T|) weakly depends on r, g cannot be a function of 77 
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alone, implying that both g and K also depend weakly on r. Taking these weak 
radial dependences of g and K into account as a perturbation, we solve equations 
(PT|) and (|^ as eigen- value problems, following Silbergleit et al. (2001). 



4 Eigen- Value Problems 

First we consider the vertical eigen- value problems. Next, the eigen- value prob- 
lems in the radial direction are considered. 



4.1 Vertical Eigen- Value Problems 

First, we solve vertical eigen- value problem, using equation (1211) . Since we are 
considering the oscillations where the vertical motions predominate over hori- 
zontal ones, we can take K to be zero in the lowest order of approximations 
(Silbergleit et al. 2001), and neglect the radial dependence oi ui^H^ /c^q, adopt- 
ing its value at a radius, say Vc- Then, equation (PT|) leads to 

(l-.^)0-2A.,| + 2iv(^)^,.O, (23) 

where the subscript c denoted the value at the radius Vc. At this stage, is 
arbitrary, but later, it is determined by an eigen-value problem in the radial 
direction. The results show that Tc can be regarded as the outer capture radius 
of the oscillations. 

Equation ([23l) is solved as an eigen-value problem with the boundary condi- 
tion g{il) < oo. The eigenfunctions are then given by the Gegenbauer polyno- 
mials, (where n = 1,2,3...), i.e., 

gin) = c^An) (24) 

and the eigenvalues are specified by 

where n is the positive integer specifying the node number of g and 

A = iV - i. (26) 
Equations (^5]) and (^51 lead to 

= vi/„l7i„ (27) 

where is given by equation ([5|), and equation (|27|) is nothing but relation ^ 
at the capture radius. The explicit form of g is 

77 for n — 1 

l{r]) (X { 1 - (1 + 2N)r]^ for n = 2 (28) 
7 - (1 + 2N/S)rf for n = 3. 
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Next, a small radial variation of Hp' is taken into account. Then, the 
equation to be solved is [see equations ([2T|) and ([23)) ] 



(l-.^)^-2iv4 + 2iV 
drj'^ drj 



5-0, (29) 



where 



e{r) = - 1. (30) 



The value of e is obviously zero at r^. . If we consider the radial dependences of 
fl± and flK, we see, after some calculation, that e > for r < Tc and e < for 
r > Tc, when uj > 0. When lo < 0, however, situation is changed and e < for 
r < Tc and e > for r > r^, as far as a» is small and |a;| is not too large. 

For the above perturbation method to be valid, the final results must guar- 
rantee that e(r) given by equation (|30p is smaller than unity in the trapped 
region. The final results really show that this is the cased Here, we briefly 
mention the case of corrugation waves (c-mode oscillations) considered by Sil- 
bergleit et al. (2001). The corrugation waves are oscillations of n = 1 and m = 1. 
In the case of n = 1, e(r) is reduced to e(r) = (D^/il^ — 1, since (w/i7j^)^ = 1 
[see equations ([25]) and p6p ]. Hence, if m = 1 we have e — {lu — il)'^ /rt]_ — 1, 
and |e(r)| ^ 1 is expected if low frequency oscillations are present. Silbergleit 
et al. (2001) showed that such oscillations are really present. This is trapping 
of c-mode oscillations. Note, however, that in this paper we are interested in 
oscillations with m ~ 2, not m = 1. 

Equation (|29p is solved by a standard perturbation method, i.e, the per- 
turbed part of g is expressed by a series of the Gegenbauer polynomials and 
thier coefficients as well as K{r) are determined from the solvability condition 
of the inhomogeneous equation, using that the Gegenbauer polynomials, C^irf), 
are orthogonal in the range of (-1,1) with the weight (1 — rf')^^^/'^. Considering 
this, we have, after some calculation, 

2NK{r) = n{n + 2A)ie(r)x(r), (31) 



where 



3r - 1 for n = 1 

X(r) = <( 2r(5r-4)/(4r2-7r + 4) for n = 2 (32) 

(3r- l)(7r-5)/(6r2 -9r + 5) for n = 3. 



4.2 Radial Eigen- Value Problem 

The results in the above subsection show that the equation to be solved as the 
eigen- value problem in the radial direction, i.e., equation (f^ . is written as 

^ For example, in the oscillations of n = 2 and rtr = (see the next subsection for the 
meaning of n^) with F = 5/3, the range of variation of e(r) in the trapped region is ~ 0.12 
(a* = 0) and ~ 0.14 (a* = 0.2). In the case of n = 3 and nr = with F = 1.3, we have 
e = 0^ 0.0135 (a, = 0) and ~ 0.0138 (a, = 0.2). 
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(Silbergleit et al. 2001), using equation ^ and (PT|) . 



UJ dr 



dhi 



ex r 



^h^^O. (33) 



2r 



Now, we introduce a new independent variable T(r) defined by 



r(r)^ l^ ^'^^'l f^^'^ dr', ^ r(r,). (34) 



a)(r') 

Then, equation (I33p is written in the form 



where 



Qhi = 0, (35) 



Equations ([55]) and (1551) show that the propagation region of oscihations is the 
region where Q > 0. The region is that of e > 0, which is found to be inside of 
rc, when w > 0. In the case of w < 0, the propagation region is outside Vc and 
oscillations are not trapped. 

Silbergleit et al. (2001) solved equation ([55]) by a standard WKB method 
with relevant boundary conditions. The WKB approximation shows that the 
solution of equation (P5|) can be represented as 

hi oc Q-i/4(t)cos [$(t) - $e] (37) 

in the whole capture region < r < Tc , except small vicinities of its boundaries 
of r = and t = t^. Here, <I>(r) is defined by 

m = r Q''\r')dr' = r QV2(/)^:Vl_^rf,'^ (38) 

and $c is a constant to be determined by boundary conditions. Concerning the 
outer boundary condition, they take into account that the capture radius, Tc, 
is a turning point of equation ([35]) since the sign of e changes there. The inner 
boundary condition adopted is vanishing of an arbitrary combination of hi and 
dhi/dr at n. In their treatment, n is taken at the marginary stable radius and 
1/cso is assumed to have a weak singularity at the radius as 1/cso oc (r — n)"'' 
with an arbitrary parameter /i. In realistic disks, however, the inner edge is 
not a singularity but the temperature and density continue smoothly inward. 
Considering this and for simplicity, we adopt /i = 0. Furthermore, we consider 
only the case of /ii = or dhi/dr = at the inner boundary. Then, their results 
of WKB analyses show that the trapping condition is 

ni/2w^ _ / + 1/4) for dhi/dr^Q . . 

^ "^^-X 7T{nr + 3/4) for = 0, 
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where nr{= 0, 1, 2, ...) is zero or a positive integer specifying the node number 
of hi in the radial direction. The constant is also determined as 

r for dhi/dr = , 
' \ tt/2 for hi=0. ^ ' 

For a given set of parameters, including spin parameter a,, polytropic index 
N , and mass of neutron stars, M, any solution of equation p9p specifies Tc, 
which gives uj of the trapped oscillation through equation (j25p . In other words, 
uj and Tc are related by equation (j25p . i.e., = a;(rc) or rc = rc{uj). Then, the 
trapping condition determines or to as functions of such parameters as a* and 
N. 



5 Numerical Results 

To obtain numerical values of the frequency, uj, and the capture radius, rc, of 
trapped oscillations, we must specify the radial distribution of acoustic speed, 
i.e., Cs(){r). Since the final results of numerical calculations depend only weakly 
on the radial dependence of Cso, we adopt the temperature distribution in the 
standard disk, where gas pressure dominates over radiation pressure and opacity 
mainly comes from the free-free processes, which is (e.g., Kato et al. 2008) 

= 1.83 X 10i6r(am)-i/^h3/5r-9/i° cm^ s-^, (41) 

where a is the conventional viscosity parameter, m(= M/Mojf| and m = 
M /Merit, Mcvit being the critical mass-flow rate defined by 

Merit = ^ = 1-40 X 10"m g s~\ (42) 

where Le is the Eddington luminosity. Parameters a and m affect on the 
frequencies of trapped oscillations only through the magnitude of Cgp. We adopt, 
throughout this paper, a = 0.1 and m = 0.3. Other parameters specifying the 
disk-star system are m(= M/Mq) and a*. We consider the cases of m = 2.0 
and a, = ^ 0.3. 

We only consider two-armed oscillations with one, two, or three node(s) in 
the vertical direction, i.e., n = 1, 2, or 30 Oscillations with more nodes in the 
vertical direction are less interesting from the view point of observability. The 
inner boundary of oscillations is taken at the radius of k = 0, i.e., at the radius 
of the marginally stable circular orbit. At the radius, we impose /ii = as the 

^ In this section and hereafter, m is often used to denote M/Mq without confusion with 
the azimuthal wavenumber m of oscillations. 

* In oscillations with odd number of n, hi{r,z) is antisymmetric with respect to the equa- 
torial plane. That is, in oscillations with n = 3, for example, h\ have one node between the 
equator and the disk surface except on the equator. See the functional form of g given in 
equation l|28| l. 
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boundary conditiorjl except in figure 5, where dhi/dr = is also considered as 
the boundary condition at rj, for comparison. The horizontal node number, rir, 
of oscillations we consider is mainly nj. = and supplementally n-^ = 1 and 2. 

Figures 1 and 2 are the propagation diagrams for oscillations of n = 1 and 2 
(figure 1) and n = 3 (figure 2), respectively. Only the oscillations of rir = and 
r = 5/3 are shown in figure 1, but three modes of oscillations, i.e., = 0, 1, 
and 2, are shown in figure 2 for the case of F = 1.25. In the cases of figure 1, the 
propagation regions of oscillations on the frequency-radius diagram are below 
the curves of 20 — rij_ and 20 — (1 + F)^/^0_l for n ~ 1 and n = 2, respectively. 
The results of numerical calculations show that the oscillations of rt = 1 and 
rir = are trapped in the radial range shown by the upper thick horizontal 
line in figure 1. The frequency uj and the capture radius Tc are, respectively, 
oj = 864Hz and Tc — 3.53rg. Outside rc, the oscillation is spatially damped. 
The radial range of trapped oscillations with n = 2 and = is shown by the 
lower thick horizontal line in figure 1. The frequency and the capture radius in 
this case are uj — 298Hz and rr = 3.69r<,. 

Trapped oscillations of n = 3 have frequencies lower than those of n = 1 and 
2, since on the propagation diagram the curve of 20 - {3ry/^n± is below those 
of 20-0_L and 20 - (1 + F)i/20i (compare fi gures 1 and 2) . It should be noted 
that the oscillation modes with rt = 3 cannot be trapped if F > 4/3, since in 
this case the propagation region is unbounded outside and the oscillations can 
propagate away infinity as shown by arrow (see the curve of 20-(3F)i/20i for 
F = 1.45 in figure 2). In figure 2, the frequency and the radial extend of trapped 
oscillations with n = 3 are shown for F = 1.25 for three modes concerning the 
radial direction; the fundamental mode (i.e., = 0) and the first two overtones 
(i.e., nr = 1 and 2). The sets of frequency and capture radius for these three 
modes of — 0, 1, and 2 are, respectively, (40.8Hz, 4.30rg), (25.8Hz, 5.84rg), 
and (17.2Hz, 7.66rg). 

Figure 3 shows the F-dependence of the capture radius Vc- As a typical case, 
the dependence is shown for oscillations with nj. = and some values of n. It is 
noted that when F is close to 4/3, the capture radius of oscillations with n = 3 
is far outside and their frequencies are low. These characteristics become more 
prominent for oscillations with > 1, although they are not shown in figure 3 
(see figure 5). 

The frequency-F relations are summarized in figure 4 for a few modes of 
oscillations in two cases of a* = and a* = 0.1. Modes of oscillations adopted 
are n = 1, 2, and 3. In all cases rir is taken to be — 0. As mentioned 
before, the oscillations with n = 3 have low frequencies. In order to examine 
characteristics of these low frequency oscillations more in detail, the frequency- 

^ If we assume that the Lagrangian variation of pressure, i.e., 5p, vanishes at r = ri, the 
condition of /ii = at r = rj will be a better approximation than dhi/dr = at r = r; 
by the following reasons. The Lagrangian variation of pressure, 5p, can be expressed as 
5p = pohi + ^rdpo/dr + ^^dpo/dz, where and are, respectively, the radial and vertical 
displacements associated with the perturbation and related to Ur and Uz by i{(jj — mQ)^r = Ur 
and — mO)^^ = u^, respectively. If we consider that equations of motion, equations I I14II 
— 11611 . and that the pressure will not change so sharply at the ri, the major term among the 
above expression for 5p will be poh\. 
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Figure 1: Frequency - radius plane (i.e., propagation diagram) showing the 
propagation region of two-armed nearly vertical oscillations. The propagation 
region of oscillation modes with n = 1 is below the curve labelled by 2f7 — U±, 
and the trapping of the oscillations with n = 1 and rir = is shown, in the 
case of r = 5/3, by the upper thick horizontal line (frequency is ~ 864Hz and 
capture radius is ^ 3.53rg). In oscillations with n — 2, the curve specifying the 
boundary of the propagation region, i.e., 217 — (1 + r)^/^r2j_, depends on F, and 
the curve for F = 5/3 is shown. The propagation region of oscillations with 
n = 2 is below this curve in the case of F = 5/3. For F = 5/3, the trapping of 
the oscillations with n = 2 and = is shown by the lower thick horizontal 
line (frequency is 298Hz and capture radius is ^ 3.69rg). The inner boundary 
condition adopted at rj is /ii = 0. This inner boundary condition is adopted 
in all cases in this paper, except for in figure 5. The central star is assumed 
to have no spin. The mass of the central star is taken to be 2Mq in all cases 
shown in figures in this paper. 

F relation in case of n = 3 is again shown in figure 5, including cases where 
other parameter values are adopted. That is, in addition to oscillations with 
Tir = 0, oscillations with n-^ = 1 and 2 are considered in figure 5. In addition, 
the cases where dhi/dr = is adopted as the inner boundary condition at rj 
are shown by thin curves. In figure 6, the frequency - spin relation is shown for 
three modes of oscillations with n = 1 , 2, and 3, where tIj. = and some values 
of F are adopted. 

6 Discussion 

In this paper we have examined characteristics of two-armed (to — 2), nearly 
vertical oscillations, assuming that the disk consists of barotropic gases with 
polytropic index F. The parameters specifying oscillation modes are n and 
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Figure 2: The same as figure 1, except that the oscillations with n — 3 are 
considered here. The propagation region of the oscillations is below the curve 
labelled by 2n - (3r)i/2f7j_, and the curve in the case of F = 1.25 is shown. 
Trapping of three modes of oscillations with tIj. = 0, 1, and 2 are shown by three 
horizontal thick lines. The sets of frequency and capture radius for these three 
oscillation modes are, in turn, (40.8Hz, 4.30rg), (25.8Hz, 5.84rg), and (17.2Hz, 
7.66rg). For the gas with F > 4/3, the propagation region is in the outer region 
of disks, which is shown by arrow, and there is no trapping. 




Figure 3: Capture radius, Tc, as functions of F for three modes of n = 1, 2, 
and 3. The radial node number rir is taken to be zero with boundary condition 
/ii = at Ti = 3rg, the spin parameter being a* = 0. In oscillations with n = 3, 
the trapping is absent for F > 4/3, when a* = 0. 
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Figure 4: Frequcncy-r relation of trapped oscillations for some values of vertical 
node number n and spin parameter a*. The oscillations with no node in the 
radial direction (rir = 0) are considered with boundary condition of /ii = at 
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Figure 5: Frequency-F relation for oscillation modes with n = 3. Effects of 
differences of radial node number rir and of boundary condition on the frequency 
of trapped oscillations are examined. Two cases of boundary conditions, hi = 
and dhi/dr = at rj, are compared for three modes of oscillations with radial 
node number rzr = 0, 1, and 2. The thick curves are for the cases where the 
inner boundary condition is taken as h\ = 0, while the thin curves are the cases 
of dhi/dr = 0. The spin parameter a* is taken to be zero. 
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Figure 6: Frequency-spin relation for modes with n = 1, 2, and 3. The radial 
node number is taken to be — 0. 



rir, where n{— 1,2,3...) is the node number of hi in the vertical direction, 
and n,:{= 0,1,2,...) is the node number of hi in the horizontal direction. An 
additional important parameter is F, which is especially important in oscillations 
with n = 3. Main results obtained are i) there are oscillation modes that are 
trapped in the inner region of disks, and ii) their frequencies depend on modes 
of oscillations (n, n^) and F, and cover a wide range of frequency. That is, the 
trapped oscillations of n = 1 and 2 have frequencies of the order of kHz QPOs, 
and those of rt = 3 arc in the frequency range of the horizontal-branch QPOs 
(HBQPOs). 

One of interesting characteristics of the oscillations is that their frequencies 
can change rather widely by change of disk structure, which is distinct from the 
g-mode oscillations trapped around the radius of K^ax (Okazaki et al. 1987), 
where Kmax is the maximum value of the epicyclic frequency. That is, as shown 
in figure 4, the frequencies of trapped oscillations decrease by increase of F, in 
the n = 2 and n = 3 modes. In the n = 1 modes, however, their frequencies 
are insensitive to a variation of F (see figure 4) . This can be easily understood 
if we consider the F-dependence of the boundary curve specifying the capture 
radius on the propagation diagram [see figures 1 and 2, and also equations © 
and 0]. 

In this paper we did not quantitatively consider the effects of Cso{r) on 
frequency. An increase of Cso without any change of other parameters leads to 
decrease of frequency of trapped oscillations. The reason is that an increase of 
Cso decreases Q. Hence, to satisfy the trapping condition ((39|) . an increase of Tc 
is necessary, which leads to decrease of frequency (see figures 1 and 2). In table 
1, the effects of changes of various parameter values on frequencies of trapped 
oscillations are summarized. 

One of important problems remained is whether the nearly vertical oscilla- 
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Table 1: table 1. 



parameters 










frequency 


parameters of oscillations 


increase 


of 


n (except for n ~ 1) 


decrease 




increase 


of 






decrease 


disk parameters 


increase 


of 


r (except for n = 1) 


decrease 




increase 


of 


CsO 




decrease 


parameters of stars 


increase 


of 


m = 


(M/Mo) 


decrease 




increase 


of 


a* 




increase 



tions considered in this paper can be really excited on disks. Two possibilities 
will be conceivable. One is the excitation by the process of viscous overstability 
of oscillations (Kato 1978). This, however, might be inefficient to some types of 
oscillations, especially to oscillations with large n. Another and more promising 
process is the stochastic excitation of oscillations by turbulence, developed first 
by Goldreich and Keely (1977a, b). This is known as the excitaion process of 
solar and stellar non-radial oscillations. This process is better than the former 
in the sense that it will be able to excite many types of oscillations, without no 
particular selection concerning the forms of eigenf unctions. 

So far, we did not discuss the effects of corotation resonance on the present 
trapped oscillations. It is known that non-axisymmetric g-mode oscillations are 
generally damped by corotation resonance (Kato 2003; Li et al. 2003; Latter 
& Balbus 2009; see also Silbergleit & Wagoner 2008). The c-mode oscillations, 
which are non-axisymmetric, are also damped by the resonance (Tsang & Lai 
2009). These results are related to the fact that these oscillations have node(s) in 
the vertical direction. The vertical oscillations considered in this paper have also 
node(s) in the vertical direction. The oscillations, however, have an important 
difference from g- and c-mode oscillations. That is, in the case of two-armed 
nearly vertical oscillations, the radius of coratation resonance, i.e., the radius 
where w = 2ri is realized, is far outside the propagation region (see the curve 
of 2fl in figure 1; the curve of 2J7 is not shown in figure 2, since it is outside 
the diagram). In the evanescent region the wave amplitude is spatially damed 
exponentially. Hence, the effects of corotation damping are negligible in the 
present oscillations. 

It will be important to note here that the effects of the general relativity 
are not essential in the trapping of the present vertical oscillations, except for 
the modes with n = 1. The propagation region of the vertical oscillations is 
specified by inequality ([8]). In deriving inequality ([8]), the fact that 5'„f25. 
larger than is adopted. In the case of oscillation modes of n = 1, the general 
relativity is necessary to guarantee this, since = 1. In oscillations with n > 2, 
however, ^n^^ > is guarantees even in Newtonian disks. Other important 
ingredients for presence of trapped oscillations are the presence of inner edge 
of disks where waves are reflected back and a monotonical decrease of angular 
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velocity of rotation, n(r), outwards. This point is different from trapping of g- 
mode oscillations, since in trapping of g-mode oscillations, the general relativity 
is essential. If we want to describe the QPOs in disks extending from black-hole 
or neutron star systems to dwarf-novae systems by a common mechanism, the 
nearly vertical oscillations considered here will be one of good candidates, since 
in the disks of dwarf- novae the effects of the general relativity are minor. 

We have shown that there are various modes of nearly vertical oscillations 
trapped in the inner region of disks. From the observational points of view, 
however, the oscillations with n — I and n = 2 may not be so interesting, since 
their trapped regions are too narrow as shown in figures 1 and 3, and may not 
be observed with large amplitudes, although the trapped regions are close to 
the inner edge of disks. Compared with them, the oscillations with n — 3 will 
be of interest, since their trapped region is wide as shown in figures 2 and 3. In 
this sense, the n = 3 modes will be one of possible candidates of low frequency 
QPOs such as horizontal- branch QPOs (HB QPOs). Further discussion on this 
direction will be worthwhile. Related to this, there are some issues to be noted 
here. 

First, we have found that low frequency oscillations with n = 3 are present 
only in disks with F < 4/3. Hence, one may imagine that such oscillations 
cannot be expected in neutron-star disks, since the standard disk model shows 
that the main part of the disks is gas-pressure-dominated. We should notice, 
however, that in disks with high accretion rate (accretion rate is close to or 
higher than the critical accretion rate defined by the Eddington luminosity), 
a radiation-pressure-dominated region appears in the innermost part of disks. 
In particular, in slim disks the whole region is radiation-pressure-dominated. 
Furthermore, we should notice that it is unnecessary to regard the index F as 
the same as the ratio of the specific heats, 7, of the disk gas. In practice, we can 
expect F < 7 by the following reasons. Due to radiation from a hot corona, disks 
have tendency to approach isothermal disks in the vertical direction. Further- 
more, in low frequency oscillations, radiative heat transport tends to make the 
oscillations isothermal, although as another effect it may dampen oscillations. 

Second, the oscillations with n = 3 are trapped in a rather wide region,, 
i.e., the radial wavelength is long. Hence, it will be necessary to take into 
account some terms neglected in deriving equation (jl9[) . including the terms 
of relativistic corrections, when we want to do more quantitative estimate of 
oscillation frequencies. 

Finally, we should mention a difference between the low frequency oscil- 
lations with n = 3 discussed in this paper and the low-frequency one-armed 
corrugation waves (Kato 1989, Silbergleit et al. 2001). The latter oscillations 
are a kind of warps or tilts, and roughly an incompressible deformation of disks, 
while the former are compressible oscillations. 

The author thanks the referee for valuable comments with careful reading 
of the original version. 
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